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Abstract 

The macroturbulence instability observed in fluxline systems during remag- 
netization of superconductors is explained. It is shown that when a region 
with flux is invaded by antiflux the interface can become unstable if there is 
a relative tangential flux motion. This condition occurs at the interface when 
the viscosity is anisotropic, e.g., due to flux guiding by twin boundaries in 
crystals. The phenomenon is similar to the instability of the tangential dis- 
continuity in classical hydrodynamics. The obtained results are supported by 
magneto-optical observations of flux distribution on the surface of a YBCO 
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single crystal with twins. 



Typeset using REVT^ 



The interest in the physics of the vortex state in type-II superconductors increased sig- 
nificantly after the discovery of high-T^ superconductivity (HTS). A main reason for the 
renewed attention is the observation of many novel nontrivial phenomena occurring in the 
vortex matter of the HTS materials. The perhaps most dramatic of these phenomena is 
the turbulence instability of the vortex-antivortex interface, which was observed in 1-2-3 
systems using magneto-optical (MO) imaging It consists in the following. When 

magnetic flux is trapped in the superconductor and a moderate fleld of reverse direction is 
subsequently applied, a boundary of zero flux density will separate regions containing flux 
and antiflux. In some temperature and fleld range this flux-antiflux distribution can display 
unstable behavior characterized by an irregular propagation of the boundary where flnger- 
like patterns develop. This contrasts strongly the regular propagation of the flux front during 
virgin fleld penetration, when only one flux polarity is present in the sample. It is clear that 
the instability, often called macroturbulence in the literature, cannot be understood within 
the frameworks of the critical state model or conventional models for flux relaxation 

An attempt to explain this remarkable behavior of the flux-antiflux interface was made 
in 1^ , where the instability was attributed to a thermal wave generated by local heat release 
in the vortex-antivortex annihilation. Unfortunately, this mechanism can hardly be accepted. 
Indeed, the vortex energy consists of two terms, where one is the magnetic part related to 
the magnetic fleld of the vortices. The other part is stored within the vortex core and 
represents the condensation energy. The magnetic energy of a vortex-antivortex pair is 
dissipated as Joule heat as they are getting close to each other, but before the annihilation 
takes place. This energy, which equals JE where J is the current density and E is the 
electric fleld generated by the vortex motion, is dissipates not only at the interface, but 
in the bulk of the sample. Hence, only the release of the core energy in the process of 
annihilation is concentrated near the interface. A simple calculation shows that the core 
energy is much smaller than the magnetic part and will cause only a negligible rise of 
the sample temperature 0], thereby ruling out that thermal effects are responsible for the 
instability. 
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In this Letter we present an explanation focussing on the experimental fact that the insta- 
bility was reported only for YBa2Cu307_5 and other 1-2-3 single crystals. These materials 
are characterized by the existence of two twin boundary systems oriented orthogonally to 
each other in the ab plane. The vortex motion in such crystals is governed by a pronounced 
guiding effect [§,0, and the viscosity coefficient for motion along the twin boundaries should 
be much smaller than across them. This anisotropy gives rise to vortex motion with a veloc- 
ity component normal to the Lorentz force. The vortices and antivortices are then forced to 
move towards each other along the interface, where the tangential component of the velocity 
becomes discontinuous. It is well known that hydrodynamical flow under such conditions 



can be unstable [jTO[, and we show here that a purely hydrodynamic theory of the vortex- 
antivortex system with anisotropic viscosity can explain the origin of the macroturbulence. 

Consider an infinite superconducting plate of thickness 2d with the external magnetic 
field H oriented parallel to the sample surface along the 2;-axis. The x-axis is perpendicular 
to the plate and x = in the center. Let H first increase and then be lowered through 
zero to a negative value. Then two kinds of vortices will exist in the sample; one with 
field direction along the positive 2;-axis {vortices) and one directed oppositely {antivortices). 
From the symmetry of the problem it is sufficient to consider only the region < x < 
and Fig. 1 shows schematically the distributions Ni{x) and N2{x) of vortex and antivortex 
densities. 

The densities Ni and N2 satisfy the continuity equation, 

C\ TV r 

^ + div(iV,V;) = 0, a = 1,2 (1) 

where Va are the vortex and antivortex velocities. In the regime of anisotropic viscous flow 
they are related to the Lorentz driving force by, 

VikNaVc^k = Flu Fl = -Bx J, (2) 

c 



B = nJo, J=^VxB. 
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Here r]ik is the symmetrical tensor of the anisotropic viscosity and $o is the magnetic flux 
quantum. Eq. can be rewritten as, 

where jik is the dimensionless tensor of the inverse viscosity, ri~f} = T'jik, and the coefficient 
r is chosen so that the principal values of 'jik are unity and e, which satisfies < £ < 1. The 
case of a strong guiding effect corresponds to e — 0. Note that F increases rapidly with the 
temperature due to thermal depinning of the vortices. 

To solve the problem one must also formulate boundary conditions at the vortex- 
antivortex interface. Generally the position of the interface, x = xq, depends on the y- 
coordinate and time t as it moves with a velocity U. The first condition is that the total 
fiux of vortices and antivortices through the interface vanishes, 

Nl{V, - U)n + N2{V2 - U)n = 0. (4) 

Secondly, both vortex polarities annihilate at the interface with the rate proportional to the 
product of their densities, 

N,{Vr-U)^ = RN^N2. (5) 

The parameter R can depend on the vortex densities and velocities but, for simplicity, we 
assume here that it is a phenomenological constant. A similar model for the annihilation 
process was used in Ref. [0. However, contrary to [0, we consider the parameter R to be 
defined by the microscopic Meissner current of individual vortices, which is much greater 
than the macroscopic currents J . Therefore, the relative velocity of annihilating vortices 
and antivortices is much greater than the hydrodynamic velocity Va- As a result, the region 
where vortices and antivortices coexist and annihilation takes place is very narrow, and can 
be represented by the surface x = xo{y,t). Finally, we take the magnetic induction to be 

A^i = N2, (6) 
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at X = xo{y,t). This condition is also a direct consequence of Eq. (D if the vortex and 
antivortex densities on the interface were equal to zero initially. 

As a first step, we determine what the model predicts for the unperturbed distribution 
profiles A^i and N2. In this case the vortex-antivortex interface is the plane x = xo(t) which 
moves with the velocity U = dxo(t)/dt. Since a simple solution with constant U ^ could 
not be found, we consider only the stationary profile with U = 0. With dNi/dt = dN2/dt = 
0, the Eqs. (|), give. 



N,{x)=N2{d) 



Xq 



d/2r — x\ 



1/2 



N2{x) = N2{d) 



d + d/2r — Xq J 
N2{d) = Ho/%, (7) 



' X + d/2r - xoY^'^ AnRd 



^d + d/2r-xoJ ' T<^hxx' 
In the following, we assume r ^ 1 since the rate of the annihilation is fast and the viscosity 

is not small. Therefore, at the interface the vortex densities are relatively small, Ni = N2 ^ 

N2{d)r^^^'^ while the spatial derivatives of Ni and N2 at x = xq are large, 

N!^ = -N[ = N^r/d, N'^ = N'; = -N^{r/dY. (8) 

Then, also the velocities Vax,y{x = xo) are relatively high, since they are proportional to 
r^/^ near the interface. 

If one cannot provide fixed values of the vortex and antivortex densities A'^i(O) and 
N2{d) at the boundaries, the discussed stationary profiles cannot be realized. In practice 
the external magnetic field defines the density of the antivortices only. Because of this 
asymmetry and the vortex-antivortex annihilation, the total number of the vortices will 
decrease with time and the interface moves towards the middle of the sample. To simplify 
the problem, we assume the interface velocity U to be much smaller than the vortex velocity 
Va, i. e. we consider the problem in the quasi- stationary regime. 

To investigate the stability of the interface with respect to small perturbations, it is 
suitable to introduce the following dimensionless variables, 
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na = Na/Na{Xo), T = t/to, to - ^^olxx 



47ri?2Ar„(xo)' 

e = x/L, C = y/L, L = ^^^ = d/r. (9) 

Normalization using the time-dependent Na{xo{t)) is allowed here since we assume that 
the instability develops much faster than noticeable changes occur in Na{xo). We seek for 
perturbations in the vortex and antivortex densities of the form 

ria = + /„exp[Ar + ikC + p^{^ - Ur))]- (10) 

The linearized boundary conditions should be written on the perturbed interface, 

e = UC, r) = Ur) + S^expiikC + Ar), (11) 

with the normal unit vector 

17= {1, -ikdaCr)). (12) 
It follows directly from Eq. (^) that 

= (/i - /2)/2. (13) 

Equations (|l]) give the expressions for the parameters pi and p2. Substituting them and 
Eqs. (j^), (|13D into Eqs. (H), (||), we obtain two linear algebraic homogeneous equations 
relating fi and f2- Demanding the determinant to vanish and omitting terms of the order of 
{Uto/Ly <^ 1, one obtains the following dispersion equation for the increment A at different 
wave numbers k, 

X = n^ -en'^ - 2isK -1-b . (14) 



Here 



k\au\ 'j^y to , . 

K = -^7—, a = — , u = U—, (15) 

2 7xx -tv 



e = Ae/{auY, s = sign(aM), 2b = —n[ — n' 
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and is a root with KeQ > of the equation 

Q'^ + + Q\-eK^ + 2 - 2b) 
-Q{2eK^ + isK + 46) - Ssin = 0. (16) 

Shown in Fig. 2 is the dependence of the increment ReA on the dimensionless wave 
number k for different e. The parameter b is set to unity as this follows from the quasi- 
stationarity condition u 0. The curves demonstrate that a positive increment exists, i.e. 
the planar interface becomes unstable, when e < ec = 0.019. The instability occurs at not 
too small values of the anisotropy and velocity of the flux-antiflux interface. Furthermore, 
it is characterised by a temporal scale given by the largest ReA, A^, which occurs at finite 
K, = K,rn- For € Sufficiently small, is much greater than unity and one finds from Eq. (p!6D 



A^ = l/(4v^) - 2, Km = \l^f2e'l\ (17) 

Thus, we have identified that an instability occurs if e < ec, which in dimensional notation 
is expressed as 



f/tan^ 



(18) 



2RNx{xq)_ 

where Q is the angle between the direction of flux guiding and the flux-antiflux interface. The 
above analysis allows us to understand why the instability of the vortex-antivortex interface 
has been observed only in crystals of the 1-2-3 system, where a pronounced guiding effect 
is expected due to the twin boundaries. 

As an illustration we show in Fig. 3 turbulent behaviour observed in an optimally doped 
YBCO single crystal containing a substantial amount of twinning, see (a). The crystal has 
a rectangular shape in the ab plane and measures 1 mm along the longest edge. Shown 
in (b) is the MO-image of flux penetration in an external fleld of 100 mT applied after 
zero-fleld-cooling to 45 K. One sees that the fleld penetrates predominantly from some large 
twin boundaries located at the bright core of the lines that make 45°-angles with the edges 
of the crystal. The dark area in the centre of the sample is the region not being penetrated 
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by a 100 mT applied field at this temperature. In (c) the temperature is raised to 67 K, 
and the MO-image was recorded in the remanent state after first applying 100 mT. The 
bright "aura" around the crystal is here the return field of the flux trapped in the central 
part. Note that this reverse field partly penetrates the sample near the edge. A distinct line 
can be seen as a dark band going around the crystal just on the inside of the edge. This 
band is the annihilation zone, which divides the crystal into two opposite magnetic domains. 
The macroturbulence is here seen as a meandering of the annihilation zone. Adjacent to 
the zone one can find small areas with increased flux density, see magnified view in (d). 
Time resolved measurements show that the zone develops in a highly dynamic manner 
where abrupt redistributions of flux often occur. By adding an external reverse field to the 
remanent state the annihilation zone is pushed further into the crystal and the dynamical 
features becomes even more spectacular. 

In this crystal the turbulent behavior was observed in the interval 25-75 K. As the 
temperature increases the dynamics of the fiux/antifiux interface becomes increasingly rapid. 
However, above 75 K there again appears to be no irregular behavior of the interface. 

The existence of the instability only in a definite temperature region finds a simple 
rationalization within our model. At low temperatures, the viscosity increases exponentially, 
and the characteristic spatial scale L, in Eq. (|), decreases correspondingly and becomes 
comparable to or less than the twin-boundary spacing. As a result, the guiding effect is 
suppressed and the instabihty disappears. On the other hand, at temperatures close to 
Tc the flux guiding is no longer effective due to thermal activation of the vortices. It is 
remarkable, and in full support of our model, that in the present heavily twinned crystal 
the turbulence occurs down to much lower temperatures than found in previous studies of 
similar crystals with only little twinning 



In conclusion, we have presented a theoretical model which gives a fully consistent de- 
scription of the macroturbulence phenomenon we observe in HTS materials under certain 
conditions. This work is supported by INTAS/RFBR and RFBR, grants IR-97-1394 and 
00-02-17145, 00-02-18032, and the Research Council of Norway. 
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FIGURES 




FIG. 1. Flux distribution in one half of an infinite slab (|a:| < d) containing trapped vortices 
of density Ni{x) in the central region \x\ < xq, and antivortices of density N2{x) penetrating from 
the outside. The other symbols are defined in the text. 
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FIG. 3. (a) - Polarized light image showing twin domains in a small area on the crystal. The 
arrow indicates the sample edge, and the scale bar is 50 /xm long, (b)-(d) are magneto-optical 
images where the brightness represents the magnitude of B's component normal to the surface, 
(b) - Applied field of Ba = 100 mT at T = 45 K. (c) - Remanent state after full flux penetration 
at T = 67 K. (d) - magnified view of the area marked in (c) 
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